Kinematics of electrons near a Van Hove singularity 
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A two dimensional electronic system, where the Fermi surface is close to a Van Hove singularity, 
shows a variety of weak coupling instabilities, and it is a convenient model to study the interplay 
between antiferromagnetism and anisotropic superconductivity. We present a detailed analysis of 
the kinematics of the electron scattering in this model. The similitudes, and differences, between a 
standard Renormalization Group approach and previous work based on parquet summations of log 2 
divergences are analyzed, with emphasis on the underlying physical processes. General properties 
of the phase diagram are discussed. 
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Since the discovery of high-T c superconductivity, the 
study of interacting electrons in two dimensions has been 
a topic of wide interest. A great variety of numerical 
and analytical techniques have been developed. Progress 
has been relatively slow, and there is little consensus on 
the phase diagram of even the simplest models. In the 
following, we revisit one of the most extensively stud- 
ied models, interacting electrons in a single band whose 
Fermi surface lies close to a Van Hove singularity. Inter- 
est in this model arose rapidly after the two dimensional 
character of the electron bands in the cuprates was estab- 
lished 0-^]. Since then, a number of different analytical 
approaches have been tried [|s| — P~0|| ■ Its interest increased 
when photoemission experiments showed the existence of 
almost dispersionless regions near the Fermi surface [ pd| , 
leading to renewed efforts to understand the mechanism 
for superconductivity using this feature. 

A few facts seem to be well established about the 
model: i) It is unstable, even at weak coupling, mak- 
ing it similar to the one dimensional Luttinger liquids 
and ii) Among the possible instabilities, antiferromag- 
netism and d-wave superconductivity seem to be the 
most likely candidates in a wide range of parameters, 
although ferromagnetism can also exist These 
features (apart from ferromagnetism) where already dis- 
cussed in the early papers on the model. It is thus sur- 
prising the little progress made in the intervening period. 

The purpose of the present paper is to disentangle the 
various factors which complicate the study of the model, 
and to study the interplay and coexistence between su- 
perconductivity and antiferromagnetism. As discussed 
in detail below, the apparent similarities with one di- 
mensional systems have lead to approximations which 
miss important features of the kinematics of electrons in 
two dimensions. This neglect, in turn, gives rise to the 
main difficulty in the study of the model, the fact that 
the perturbative divergences mentioned before seem to go 
like the square of the logarithm of the cutoff. This fea- 
ture, never found in a renormalizable field theory in the 
usual sense, has been dealt with by parquet summation 
techniques, using the analogy with one dimensional sys- 



tems. While these approximations have their own merits, 
they cannot be considered a proper implementation of a 
Renormalization Group procedure. In the following, we 
extend previous results, and present a detailed analysis 
of the renormalization of interactions in this model, us- 
ing the by now standard approach presented in [^4|-^6[. 
The present work shows that the Van Hove singularity 
problem can be treated by techniques which are easily 
extended to generic two dimensional problems, allowing 
us to extrapolate the phase diagram to different ranges 
of parameters and fillings. 
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FIG. 1. Bare couplings. Full and dashed lines stand for 
propagators at the two inequivalent saddle points. 

We focus on the weak coupling instabilities of the 
model. As mentioned earlier, a perturbation analysis 
shows that, unlike models with isotropic Fermi surfaces, 
the RG flow of the electron-electron interactions is not 
trivial. We consider only the electronic states lying near 
the two inequivalent saddle points of a square lattice. Fol- 
lowing preceding analyses, we classify the couplings into 
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four groups, shown in Fig. |], denoted as intrasingular- 
ity, intersingularity, backwards and Umklapp scattering. 
In addition to this, the conservation of total momentum, 
and the requirement that the scattered electrons lie close 
to the Fermi surface induce new constraints. The pair 
of outgoing momenta is uniquely determined by the pair 
of incoming momenta, except when they add up to zero 
jl4],[l6| . This allows us to define the Cooper pair channel 
in this latter instance, besides the forward and exchange 
channels, that open up when the momentum transfer to 
any of the outgoing particles vanishes (modulo (n, ir) in 
the present model). 
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FIG. 2. Diagrams which correct the bare couplings at one 

loop. 

The lowest order corrections to the bare couplings are 
shown in Fig. |^. In a conventional Fermi liquid, only the 
diagram (a) in Fig. || depends logarithmically on the cut- 
off, leading to the BCS instability. In the present case all 
diagrams show logarithmic dependences, if the momen- 
tum of the excitations which are integrated over is close 
to zero. The ladder diagram (a) in Fig. || shows a log 2 (A) 
dependence, while the others behave as log(A). We are 
not considering the possibility of perfect nesting, when 
the Fermi surfaces at the two singularities are parallel 

0- 

We first discuss the diagrams (b), (c) and (d) in Fig. 
0, The restriction on the intermediate momenta implies 
that the divergences appear when k\ m k$. Thus, di- 
agrams (b) and (c) give non trivial corrections in the 
forward channel, and diagram (d) modifies the exchange 
channel. Note however that, when all spins are parallel, 
the forward and exchange channels can only be distin- 
guished by the momentum transfer involved, as they lead 
to the same final state. Physical observables depend only 
on the difference between them, like the direct and 2kp 
scattering between electrons with parallel spin in one di- 
mension. In the following, we make use of this fact and 
redefine the parallel forward channel as this combination. 
In this way, the exchange channel only has antiparal- 
lel couplings. Furthermore, the forward channel is only 
renormalized by couplings in the forward channel. The 
same applies to the exchange and Cooper pair channels. 

With the mentioned redefinition of the forward parallel 
couplings, the net effect of the renormalization is given by 
diagram (b), which tends to screen the bare interaction. 
The RG equations in the forward channel are: 
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where c, c' are the prefactors of the polarizabilities at zero 
and Q = {%, w) momentum transfer. A formally equiv- 
alent equation is obtained for the couplings in the ex- 
change channel, except that, in this case, the signs in the 
r. h. s. of Eq. (^) are positive, and repulsive couplings 
increase upon scaling. The equations in ([!]) can be fur- 
ther simplified. For instance, we can add and substract 
the upper four equations, to obtain: 
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FIG. 3. Generic flow of the couplings. The flow direction is 
towards the left in the forward channel, for the combination 
of couplings in Eq. (j^). A similar flow diagram can be drawn 
for the Cooper pair channel, with the replacement of un , iii 
by the Cooper pair couplings ucintr a ± , wcumkx • 

These equations can be integrated analytically. The 
main features of the flow are depicted in Fig. |[ The 
flow for the forward channel is towards negative values 
of the abscissas, the opposite as for the exchange channel. 
The hamiltonian flows towards different fixed points for 
different initial conditions, defining different universality 
classes. When the initial couplings are spin independent 
P,p|pO|] , the flow starts at the diagonal. The couplings in 
the forward channel vanish in the infrared. If the bare 
hamiltonian is the Hubbard model in its usual form, the 
initial parallel couplings are zero. Then, the couplings 
must flow towards a strong coupling regime. If only the 
forward and exchange channels are considered, the lead- 
ing instability is towards a broken symmetry state with 
antiferro- or ferromagnetic order, depending on whether 
c' or c is larger. It is easy to show that the (attrac- 
tive) divergent flow of the parallel couplings in the for- 
ward channel matches the repulsive flow of the antiparal- 
lel couplings in the exchange channel. The two channels 
renormalize the parallel and transverse spin-spin correla- 
tions, respectively, so that the flow maintains the original 
SU(2) symmetry of the microscopic model. In the follow- 
ing, we consider the modifications to this picture induced 
by the inclusion of the Cooper channel, which has to be 
considered independently. 

The main source of difficulties, the log 2 divergences in 
the perturbation expansion, is found in the Cooper chan- 
nel. A number of previous works have dealt with them 
by scaling the couplings as function of log 2 (A), instead of 
the more standard log(A). Let us consider in some detail 
the implications of this procedure. In a renormalizablc 
theory in the usual sense, an effective hamiltonian can be 
defined at each scale, A, with dimensionless interactions, 
{v.i}, that depend on A. The elegance of the method re- 
sides in the fact that there is no explicit dependence on 
A elsewhere. The most general equation for the flow of 
the couplings is: 

where the functions / are dimensionless and A stands 
for additional low energy scales. The statement that the 
model is renormalizable implies that the limit A/ A — » 
can be safely taken. On the other hand, a scaling with 
log 2 implies that the r. h. s. of Eq. (||) is of the type: 

/(W'£) =/({M)l°g(£) (4) 

The dependence on the low energy scales represented by 
A cannot be avoided. A scaling in terms of log 2 is not 
equivalent to a conventional RG calculation, as it ignores 



the low energy scales present in the flow equations. It 
resembles more a parquet summation of leading diver- 
gences. In one dimension, where only logarithmic diver- 
gences appear, both approaches are equivalent. This is 
not the case in the present problem. It can be shown 
H that a log 2 divergence, when inserted into a self en- 
ergy diagram, leads to further log 2 divergences in the 
one electron propagator. The usual wavefunction renor- 
malization cannot take care of this divergence, signalling 
again the existence of hidden scales and/or new, non local 
interactions. 

In order to avoid the log 2 divergences in the Cooper 
channel, we move the chemical potential slightly away 
from the singularity. It is easy to show that the diver- 
gence in the Cooper channel becomes log(A/a>) log(A//x), 
where /i is the chemical potential. If [i -C A, the other 
channels are unaffected. The factor log(A/^t) is simply 
the density of states at the Fermi level, which, in prin- 
ciple, can also be scale dependent. A full treatment of 
this term requires the analysis of the flow of the chemical 
potential itself. 

We will assume that the system is in contact with an 
external reservoir which fixes the physical, fully renor- 
malized chemical potential. The effective chemical po- 
tential which must be inserted into the low energy hamil- 
tonian at a scale A is renormalizcd by a Hartree diagram 
In the present case, this diagram has a logarith- 
mic singularity, besides the usual linear dependence on 
A dictated by dimensional considerations . As a re- 
sult, one finds for the dimensionless chemical potential, 
ft = n/A: 

= ft - (uFintra + Winter) log ( — ^ J (5) 

olog(A) Vi-AV 

where a summation over parallel and perpendicular spin 
orientations is assumed, and the couplings have been 
written in a dimensionless form. This equation has a 
stable fixed point, ft*. This fixed point is attained if 
ft > 0, that is, if the chemical potential lies above the 
singularity. The existence of this fixed point is crucial 
in simplifying the treatment of the scaling in the Cooper 
channel. For sufficiently small initial couplings or at suf- 
ficiently low energies, we can replace the log(A//x) in the 
scaling of the couplings in the Cooper channel by the 
constant log (1 /ft*). Note that the previous analysis does 
not contradict the conservation of the number of elec- 
trons implied in Luttinger's theorem, as we have chosen 
to work at fixed external chemical potential. We find 
that this is the natural choice in the problem at hand. 
The restriction of the study to small patches near the 
saddle points is justified as these regions give rise to the 
main divergences in perturbation theory. However, the 
remaining sections of the Fermi surface are, at least, a 
reservoir of electrons, always in contact with the system 
studied here. 

The previous discussion shows that one of the two 
log(A) factors which appear in the Cooper channel should 
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not be included in the RG equations, provided that the 
chemical potential is shifted slightly away from the sin- 
gularity. The resulting equations are analogous to those 
in Eq. ([!]), except that the coefficients include a con- 
stant equal to log(l//i*) at the beginning of the scaling. 
Note, finally, that the existence of a finite chemical po- 
tential implies that the flow in Eqs. ([!]) ceases to be 
valid at u) ~ fi. On the other hand, the couplings reach 
a strong coupling regime at scales E c ~ Ae -1 /" where 
i*o stands for the initial dimcnsionlcss couplings. Hence, 
the phase diagram, which is determined by the couplings 
which reach first the strong coupling limit, should be in- 
dependent of fj,, provided that /i -C E c . 

We now analyze the phase diagram by identify- 
ing the couplings which diverge at the highest scale. 
A simple inspection shows that ferromagnetism pre- 
vails when c > c', associated with the divergences 

in UFintra||,"Fmtra±,UFintcr|| and UFintcr J_ ■ When 

c < c! the leading divergences arc either in the set 
Whacky, WFback_L,WFumk||, «fumki } or in the couplings 
for the Cooper channel ucintra± and M Cumk_L- The first 
set gives rise to antiferromagnetism, while the second pair 
leads to d-wave superconductivity, even if the initial cou- 
plings are repulsive, provided that uc"umk± > ucmtra± > 
as shown in Fig. ||. In the Hubbard model, these 
couplings are initially equal. Finite corrections, not con- 
sidered here, influence differently ucumk± and ucintrai, 
leading to the conditions for d-wave superconductivity. 
These finite corrections are those considered by Kohn 
and Luttinger in their pioneering work on anisotropic su- 
perconductivity from repulsive interactions (ll| . 

Thus, the phases which arise in a weak coupling treat- 
ment of the Hubbard model near a Van Hove singular- 
ity are ferromagnetism antiferromagnetism and d-wave 
superconductivity. The competition between antiferro- 
magnetism and superconductivity depends on the value 
of d and log(/i*). As both types of couplings diverge, the 
transition between the two is discontinuous ( the phase 
diagram in the strong coupling limit to which the sys- 
tem flows should be well described within a mean field 
approximation with the appropiate couplings Jl9|l). Fi- 
nally, the charge compressibility depends on the density 
of states at the Fermi level divided by one plus the sum 
of the parallel and antiparallel couplings in the intrasin- 
gularity forward channel. As discussed earlier, parallel 
couplings flow towards divergent negative values, while 
antiparallel couplings are divergent and positive. The 
sum tends to zero. Thus, the compressibility approaches 
the density of states, which, by definition, is large. We 
find no evidence of a vanishing compressibility |10| . On 
the other hand, we cannot rule out that the cancellation 
found here does not exist when the Fermi level is suf- 
ficiently far from the singularity, leading to an infinite 
compressibility and phase separation pcfl . 

In conclusion, we have analyzed in detail the kinemat- 
ics of electrons near a two dimensional Van Hove singu- 
larity. While the model has been studied extensively, the 
classification of the different allowed channels has never 



been presented. This framework, which makes connex- 
ion with standard applications of the Renormalization 
Group program to two and three dimensional electronic 
systems |Q, clarifies the different physical processes in- 
volved, and the main features of the phase diagram. The 
results can differ substantially from those obtained by 
summing, in an uncontrolled way, log 2 (A) divergences. 
Finally, it can be easily extended to address other types 
of anisotropic Fermi surfaces plLf22| . 
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